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Abstract 

We show that the size distributions of fragments created by high 
energy nuclear colhsions are remarkably well reproduced within the 
framework of a parameter free percolation model. We discuss two 
possible scenarios to explain this agreement and suggest that perco- 
lation could be an universal mechanism to explain the fragmentation 
of simple fluids. 

PACS. 36.40.Ei Phase transitions in clusters - 36.40.Qv Stability and 
fragmentation of clusters - 25.70.Mn Nuclear fragmentation 

1 Introduction 

Fragmentation phenomena concern a wide diversity of objects in nature, at 
many scales of distance and time. A natural question is to ask what is generic 
and what is specific in these phenomena. Most theoretical efforts made to 
understand fragmentation apply only to specific objects or, on the contrary, 
concern simple models with few links with reality . Moreover, experimental 
data are rather sparse and often the experimental conditions are ill 

defined. As a consequence, the question of the possible existence of universal 
fragmentation mechanisms remains an open problem f^. 
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The arguments in favour of the existence of universal mechanisms are 
of various orders. For instance, in aggregation phenomena, which can be 
considered as the opposite of fragmentation, it is possible to define univer- 
sal classes (Diffusion Limited Aggregation, Clustering of clusters...) in 
terms of the initial conditions (number of seeds) and the motion of the aggre- 
gating particles (ballistic, Brownian). The fractal structure (dimension) of 
the aggregation cluster is the fingerprint of these universal classes. In many 
fragmentation processes the experimental observation, over many orders of 
magnitude, of power law (scale invariant) fragment size distributions is an- 
other possible indication of universal classes. In this case, the value of the 
power law exponent would be the corresponding fingerprint 

Collision induced fragmentation of small fluid drops (atomic nuclei 0, 
atomic aggregates |jlO|,|ll[, liquid droplets [|12|) is a field of experimental re- 
search particularly active because it offers the best possibilities of complete 
identification of the fragmentation products. We present in this paper an 
analysis of the fragmentation of atomic nuclei in high energy collisions. We 
show that random percolation theory accounts for experimental data with- 
out any adjustable parameter and we discuss two possible explanations for 
this agreement, depending if one assumes or not that thermal equilibrium is 
reached before fragmentation occurs. It is then suggested that this percola- 
tion type fragmentation mechanism could be universal for simple fluids, i. e. 
fluids made of structureless particles interacting with short range potentials. 

The structure of the paper is the following. In section 2, a brief description 
of the experiment precedes the analysis of data. Section 3 is devoted to 
the interpretation of the results. Some final remarks and conclusions are in 
section 4. 



2 Analysis of fragmentation data 
2.1 The experiment 

The experiment was performed by the Aladin collaboration at Gesellschaft 



fiir Schwerionenschung (G.S.I.), Darmstadt [|T3[. Beams of gold ions (Au, 
Z = 79) at 600 and 1000 MeV/nucleon incident energies were used to bom- 
bard targets made of thin copper foils. Data for more than 3.10^ events were 
collected. At these high energie^, the commonly admitted scenario of the 
collision is the following: The part of the projectile that does not geometri- 
cally overlap with the target at the point of closest approach is thought to 



^High energies compared with the 8 MeV/nucleon binding energies or the 35 
McV/nucleon Fermi energies of nuclei. 
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Figure 1: Experimental fragment size distributions n{z) at four values of the 
control parameter Zbound (see text). The lines correspond to the percolation 
calculation. A log-log representation has been chosen to emphasize the power 
law behaviour at Zbound — 48. 
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Figure 2: Tesi o/ i/ie scaling assumption (Eq. 2), that requires that all the 
ratios n{z)/nc{z) follow the same curve. 



decouple from the rest and form a sub-system called the projectile specta- 
tor (PS). The size of this system and its excitation energy E* are therefore 
dependent on the impact parameter. The overlapping part of the system is 
completely vaporized and therefore does not contribute to the formation of 
fragments with nuclear charge [z) greater than one. 

In the present experiment, the Aladin device detected with very high effi- 
ciency all spectator fragments (those resulting from the decay of the unstable 
PS) with nuclear charge z > 1. Neither hydrogen isotopes nor neutrons were 
detected, event by event, with any significant efficiency. The initial size of 
the spectator system (Zps) is not experimentally measured and can only be 
inferred from the comparison with a nuclear reaction model. Empirically, 
this size can be estimated, on average, from the following relation [p!4|: 



< Zps >=25 + Z 

bound 0.004ZL„, (1) 

where Zbound is equal to the sum of all the charges of products with z > 
2. In all the studies of these experiments, the parameter Zbound is used as 
control parameter. Notice that Zhound, which is closely related to the impact 
parameter, decreases when the violence of the collision increases. 



2.2 Analysis of data 

To analyze the Aladin data, the following very simple procedure has been 
used. For a given value of Zbound, the size of the PS was deduced from eq.|l] and 
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Figure 3: Mean value of Zmax as function of Zbound for percolation cal- 
culation and experimental data from ref. |7^/ (left). Fluctuations of Zmax 
as function of Zhound for percolation calculation and experimental data from 
ref. ^ (right). 



percolation calculations on a lattice of corresponding size were performed by 
varying the bond breaking parameter (random-bond percolation). Only those 
events with the proper Zbound were kept and the fragment size distribution 
of these compared to the corresponding experimental data. In this manner 
the comparison is parameter free. Notice that by fixing a value of Zbound, the 
maximum value Zmax of z is also constrained by this value. 

Figure |l] shows, for Zbound values of 20, 48, 60 and 70, the experimental 
fragment size distributions. One observes three different regimes of fragmen- 
tation. For large Zbound, large impact parameter and low excitation energy, 
only a heavy residue and light particles are produced. This is the "evapora- 
tion" regime. For small Zbound (most violent collisions) only small fragments 
are produced and no heavy residue is left. The distribution of this "vapor- 
ization" regime decreases exponentially. Remarkably, at Zbound around 48 an 
intermediate regime exists, for which the distribution can be fitted by a power 
law n{z) ~ The mean excitation energies per particle can be estimated 

for the fragmenting systems as a function of Zbound HI] : they correspond 



to about < E* >~ 14,6,3,1 MeV per particle for Zbound = 20,48,60,70 
respectively. 

The full lines in figure 1 represent the corresponding distributions ob- 
tained from the percolation calculation. The calculation reproduces quanti- 
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Figure 4: The distribution P{zmax) of the largest fragment at four values of 
Zbound- The circles correspond to the data from ref. / (7^/ and the lines to the 
percolation calculation. 
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tatively, and over several orders of magnitude, the data in the three regime^. 
We are naturally led to associate the experimental power law behaviour at 
Zbound — 48 to the percolation critical behaviour. 

Close to the critical regime, percolation theory predicts typical scaling 
properties of the cluster size distributions |jl5[. The size distributions can be 
reduced to an universal function, /(§ ), by the following formula : 

n{z) = n,{z)f{-) (2) 

z 

where nc{z) is the size distribution at the critical point, taken here as the size 
distribution observed for Zbound = 48. The "characteristic size" z is defined 
by: 

Z = 712^/1712 (3) 

with rrik = Yl,z>2^^^'{^)i where n'{z) is the mean fragment size distribution 
obtained excluding, event by event, the largest fragment. One observes in 
figure 2 that, within the scattering of the data, this rule is rather well satisfied, 
keeping in mind that this corresponds to the ratio of two quantities varying 
over more than three orders of magnitude. 

The size of the largest fragment plays in percolation theory the role of the 
order parameter. For an infinite system, it is of infinite size in the percolating 
phase and finite in the non-percolating one. In a finite system, this transition 
is smooth, as illustrated in figure ^ This figure on the left compares, again as 
a function of Z^^oundi the experimental measured size of the heaviest fragment 
Zmax to the one obtained from the percolation calculations. The figure on 
the right shows the fluctuations of Zmax, 

2 _ ^ ^max > ~ < Zmax > / 

zmax ^ r ^ ' 

' Zmax ^ 



As expected |T^, these functions show a maximum around the "critical'' 



value of Zhound- The agreement between experiment and theory is good, 
except for the highest values Zhound where the fluctuations should vanish for 
z = 79. Very probably, either contributions of the target have been included 
in the data set or miss-identification of fragments is present. Naturally, the 
percolation model used here cannot account for either of these features. In 
figure ^ one can see that even the full distribution P{zmax) is well accounted 
for by the calculation. 



close inspection of the data for Zhound — 70 shows an excess of fragments produced 
around 2; = 30 due to the presence of fission events. This is understandably beyond the 
scope of percolation theory. 
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By any standard, the agreement observed in figures [1| to 0, between the 
experiment and the calculations is very good. In the analysis presented here, 
as already stressed, no adjustable parameters are used. 



3 Interpretation of the results 

The choice of percolation theory to analyze the experimental results is moti- 
vated by the following reasoning. A fully microscopic description of nuclear 
fragmentation is, a priori, out of scope of theory. Atomic nuclei behave in 
their ground state as small drops of Fermi liquids composed by particles 
strongly interacting mainly with a two-body, short range force. This inter- 
action is ill defined at short distances and the techniques used to solve this 
complicated many body problem are not fully under control. The description 
of collisions, using transport equations for example, is even more difficult. In 
view of the above remarks and with the stated goal of understanding the 
universal features of the data, we believe that a more fruitful point of view 
is to tackle the problem with the minimum number of assumptions. 

The simplest hypothesis would be to assume the equiprobability of all 
partitions of the integer number Z = 79. This is equivalent to a "maximum 
entropy principle" []T^|T^. However, the resulting n{z) are always exponen- 



tially decaying functions, in contradiction with experiments. 

A step further is to consider topological constraints, by considering frag- 
ments in a 3-dimensional space. Among the infinity of models one can imag- 
ine to make fragments, random-bond percolation seems a good candidate 
because of its simplicity and because it retains only the essential constraints. 
For example, the shape of the fragments is not assumed a priori. It turns out 
that, as shown in section 2, it suffices to reproduce very well the experimental 
data. How does one understand this agreement ? One can imagine at least 
two scenarios. 

In the first scenario, one idealizes the nuclear fiuid as an ensemble of par- 
ticles connected by bonds. The bond between a pair of particles is active as 
long as the magnitude of their potential energy is greater than the relative 
kinetic energy. During the collision some of these bonds break because of the 
change in the position and/or the velocity of the particles. The simplest as- 
sumption is that bonds are broken randomly, which corresponds to the usual 
uncorrelated bond percolation model |T5[. Particles connected by unbroken 
bonds form the fragments. These fragments separate rapidly, pushed away 
by the long range Coulomb force between protons. 

In a second scenario, one assumes that after the collision phase, equi- 
librium thermodynamics applies. The system expands until a "freezout" 



8 



density is reached, at which the fragments cease to interact by the strong 
nuclear attractive force and their size distribution is "frozen" . The Coulomb 
force accelerates the fragments, as before. In order to calculate its distribu- 
tion, we have to define first what a "fragment" is. In the present context, it 
seems natural to call "fragment" a self-bound ensemble of particles |19, 20[]P| . 



With these definitions, one finds in the p — T diagram a percolation line 



sometimes called the Kertesz line [^) that separates a percolating and a 



non-percolating phases. The line joins the thermodynamical critical point to 
the random bond percolation critical point. On it the fragment size distribu- 
tion is a power law n{z) ~ z"'^ , with r ~ 2.2. (see fig. 1 and 2 of ref. |^). In 
a small system like an atomic nucleus, rather than a sharp critical line, one 
finds a "critical zone" on which the n{z) approaches this behaviour (see figure 



4 of ref. [^). The similarity of these calculated n{z) with the experimental 
ones shown in figure 1 is very obvious. Therefore, by inspection of the n{z) 
alone, it is not possible to disentangle between these two scenarios. However, 
one could hope that extra information on, for example the fragment kinetic 
energies, could indicate if thermalization is present or not. 

More generally, the main difficulty in analyzing nuclear fragmentation 
data is due to the very small size of the system, which fundamentally limits 
the extraction of the universal critical exponents. Indeed, a proper charac- 
terization of the physical process requires, apart from the r exponent, the de- 



termination of the other critical exponents [|1^] associated with the moments 
TUk of the fragment size distribution. Such measurements are, in principle, 
possible for simple fluid systems (i.e. made of structureless particles subject 
to short range forces) of larger sizes such as atomic aggregates or macroscopic 
pieces of matter, such as liquid drops [1^. 

We consider as very plausible the possibility to observe this percolation 
type fragmentation in simple fluids. Indeed, the arguments that we have 
developed to explain the success of percolation theory should apply without 
restriction to any system of structureless particles interacting with a short 
range potential. In fact, attempts have been made to show experimentally 
this behaviour in the fragmentation of hydrogen aggregates In these 



experiments, the multiplicity of fragments mo is used as control parameter. 
As a function of mo the mass of the largest fragment and its fluctuation. 



evolve qualitatively as expected in percolation theory. However, in ref. |Tl 
the data are compared with a percolation system of improper size and no 
definite conclusions can be drawn. 



■^Another possibility is to impose the stability by particle evaporation only . At the 
present level of accuracy, both definitions are operationally very similar. Definition pl| is 
equivalent to the ones proposed by Hill [E3] and by Coniglio and Klein p3] 
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On the theoretical side, we are performing large scale classical molecular 
dynamics simulations of a Lennard- Jones fluid, for both the sudden disassem- 
bly of an equilibrated system and for the collisions of drops. We clearly find 
in these calculations the fragment size distributions predicted by percolation 
theory p6| . 



4 Final remarks 

Other experiments of nuclear fragmentation at high bombarding energies 
have been successfully interpreted with percolation theory |]27|- p0|j . However, 



at lower bombarding energies (30-50 MeV/nucleon) the agreement is less 
satisfactory [pT| , |3^ . The shape of the n{z) evolves qualitatively as in figure 1, 



but a closer examination shows a systematic over-production of intermediate 
size fragments. The origin of these discrepancies is still unclear. Different 
explanations can be considered: 

a) the identification of the source of the fragments is more difficult at lower 
coUisionnal energies where the separation between the projectile, the target, 
or the fused system is not as clear as in the case of high energy collisions, 
leading thus to possible contaminations between these different sub-systems. 

b) different fragmentation mechanisms could result from the smaller rela- 
tive velocities between projectile and target, particularly in the case of fusion 
at small impact parameter, inducing possible compression effects. 

These discrepancies, more generally, could just signal the fact that the 
assumption that fragmenting nuclei behave as "simple fluids" breaks down 
at these lower incident energies. 

Nuclear fragmentation experiments are often analyzed with the so-called 
Statistical Multifragmentation Models [^,0]. In brief, these models deal 



with the equilibrium thermodynamics of ensembles of spherical drops of nu- 
clear matter confined in a "freezout" volume. Drops interact with each other 
only by the long range Coulomb force and their internal partition function 
is taken from empirical mass formulas or from experiment. These models^ 
are successful in describing the above mentioned low bombarding energy ex- 
periments when fixing the size, the density and the excitation energy of the 
fragmenting source. For the high energy Aladin data similar analysis have 
been performed but with a larger set of input parameters ||28[ . 



The use of percolation theory concepts is however more comprehensive 
while much easier to handle. It provides an excellent agreement with the 
data without requiring any adjustable parameter. 



These models have been also implemented to study the fragmentation of atomic clus- 
ters isf. 
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We hope that the present results will encourage both the theoretical and 
the experimental studies of the fragmentation of simple fluids. The fragmen- 
tation by collisions of very large aggregates seems particularly promising, be- 
cause it combines the experimental possibility to detect fragments together 
with a reduction of finite size and surface corrections effects. 

We thank the Aladin collaboration at G.S.I, for allowing us to use their 
experimental data. 
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